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PROJECT  OVERVIEW 


OBJECTIVES 

The  general  objective  of  the  project  is  fundamental  mathematical  modelling  of  a  complex 
TSA  system  with  electrothermal  desorption  step,  and  with  integrated  condenser  of  the 
desorbed  compound,  which  would  enable  better  fundamental  insight  in  to  the  process  and  its 
optimization.  The  objective  of  this  stage  of  the  project  is  modeling  of  the  central  unit  of  the 
TSA  system:  the  adsorption  bed.  The  adsorbers  in  the  TSA  system  under  consideration 
consist  of  cartridge-type,  radial  flow  fixed  beds,  made  of  spirally  coiled  activated  carbon  fiber 
clot  (ACFC).  Accordingly,  the  subject  of  the  mathematical  modeling  in  this  project  is  a  single 
annular,  radial-flow,  cartridge-type  ACFC  adsorption  bed.  Models  of  different  complexity  are 
to  be  postulated  and  used  for  simulation  of  adsorption,  electroresistive  heating,  electrothermal 
desorption  and  consecutive  adsorption-desorption  in  the  described  adsorber. 


STATUS  OF  EFFORT 

Two  mathematical  model  of  the  system  under  consideration  have  been  postulated:  the  first 
assuming  uniform  adsorbent  density  throughout  the  adsorbent  bed,  and  the  second  taking  into 
account  the  structure  of  the  bed  made  of  layers  of  ACFC.  Both  models  were  obtained  as 
complex  sets  of  coupled  nonlinear  PDEs,  ODEs  and  algebraic  equations,  defining  the 
equilibrium  and  transport  parameters  of  the  system.  Numerical  procedures  for  solving  the 
model  equations  have  been  established  and  used  for  simulation  of  adsorption,  electroresistive 
heating,  electrothermal  desorption  and  consecutive  adsorption-desorption.  Nevertheless, 
comparison  with  experiments  couldn’t  be  performed  owing  to  the  lack  of  detailed 
experimental  data  corresponding  to  a  single  cartridge  adsorber.  Also,  lack  of  reliable  relations 
for  calculation  of  transport  parameters  for  an  adsorbent  bed  with  fibrous  structure  was 
observed. 


ACCOMPLISHMENTS  /NEW  FINDINGS 

Mathematical  models  for  an  annular,  radial-flow  adsorption  bed  have  been  developed  and  the 
procedures  for  their  numerical  solutions  established.  These  models  can  be  used,  not  only  for 
the  TSA  system  with  electrothermal  desorption  step,  but,  after  minor  adjustments,  for  other 
cases,  with  external  or  internal  heating  of  the  adsorption  bed,  as  well.  The  models  developed 
in  this  project  can  be  used  as  the  basic  block  for  modeling  of  the  entire  TSA  system. 


PERSONNEL  SUPPORTED 

Danijela  Antov,  MS  student 

Dr.  Menka  Petkovska,  Assoc.  Professor 
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PUBLICATIONS 

No  publications  yet. 

Two  journal  papers  in  preparation: 

A  review  paper  on  electrothermal  desorption 

A  scientific  paper  based  on  the  results  presented  in  this  report 

One  paper  accepted  for  presentation  at  the  Fundamentals  of  Adsorption  -  FOA8  Conference 
in  Arizona,  23-27  May  2004: 

-  “Electrothermal  desorption  in  an  annular  -  radial  flow  -  ACFC  adsorber  - 
Mathematical  modeling”,  by  M.  Petkovska,  D.  Antov  and  P.  Sullivan 


NEW  DISCOVERYS,  INVENTIONS  AND  PATENTS: 

None 
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DETAILED  REPORT 

1.  Introduction 

The  idea  about  regeneration  of  used  adsorbent  beds  by  direct  heating  by  passing 
electric  current  through  them  (owing  to  Joule  effect),  was  first  published  in  the  1970s  [1]. 
Desorption  process  based  on  this  principle  was  later  named  electrothermal  desorption  [2],  It 
was  recognized  as  a  perspective  way  to  perform  the  desorption  steps  of  TSA  cycles  [2],  At  the 
same  time,  fibrous  activated  carbon  was  recognized  as  a  very  convenient  adsorbent  form  for 
its  realization  [2],  In  the  last  couple  of  years  application  of  electrothermal  desorption  has 
reached  the  industrial  scale  [11,  12], 

Efflu  ent 


Condensate 


Figure  1.  Overall  schematic  of  the  ACFC  adsorption  -  rapid  electrothermal  desorption  system 

(from  Ref.  [3]) 


Thermocouple  | 


Figure  2.  Structure  of  the  adsorbers 
(from  Ref  [3]) 


This  project  is  devoted  to  modeling  of  a  novel  TSA  system 
with  electrothermal  desorption  step,  described  in  Ref.  [3],  and 
shown  in  Figure  1.  Each  adsorber  of  this  system  consists  of  two 
annular  cartridge-type  fixed  beds,  made  of  activated  carbon  fiber 
cloth  (ACFC)  (Figure  2). 

There  are  several  reasons  for  developing  fundamental 
mathematical  models  of  the  investigated  process.  This  modeling 
should  enable: 

•  Better  phenomenological  understanding  of  the  process  and 
recognizing  of  the  most  important  mechanisms  and 
phenomena; 

•  Better  insight  into  the  process  which  can  help  to  detect  it’s 
drawbacks  and  to  propose  the  necessary  improvements, 
regarding  the  necessity  for  new  measurements,  potential 
change  of  the  operational  conditions  or  even  potential 
change  of  the  design; 

•  The  possibility  for  prediction  of  the  system  behavior  at 
different  operational  conditions,  different  scale,  or  with  a 
changed  design,  which  could  substantially  lower  the  need 
for  experimental  trials; 
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•  Optimization  of  the  system  based  on  the  computational  analysis. 

Before  making  an  attempt  to  derive  a  complete  mathematical  model  of  this  complex 
system,  mathematical  modeling  of  a  single  adsorption  bed,  would  be  very  useful.  This  report 
presents  the  results  of  such  modeling. 


2.  Description  of  the  system  -  annular,  radial-flow,  cartridge-type  adsorber 

The  investigation  during  the  first  year  of  the  project  "Rigorous  mathematical 
modeling  of  an  adsorption  system  with  electrothermal  desorption  of  the  used  adsorbent",  was 
focused  on  mathematical  modeling  of  a  single-cartridge  adsorber,  for  three  processes  that 
could  be  performed  in  it:  adsorption,  electroresistive  heating  and  electrothermal  desorption  of 
a  previously  adsorbed  sorbate. 

The  annular,  radial-flow,  cartridge-type  fixed-bed  is 
schematically  shown  in  Figure  3.  It  is  formed  as  a 
cylindrical  roll  of  activated  carbon  fiber  cloth  (ACFC) 
spirally  coiled  around  a  porous  central  pipe,  used  for 
introducing  the  gas  stream  (the  stream  that  has  to  be 
purified,  during  adsorption,  or  an  inert  gas  stream,  during 
desorption).  The  gas  flows  in  the  radial  direction  through 
the  adsorber,  leaving  the  adsorber  through  an  annular  space 
around  the  adsoption  bed.  During  the  desorption  step, 
electric  current  is  passed  through  the  activated  carbon  cloth 
in  the  axial  direction,  causing  heat  generation  inside  the 
carbon  fibers  (Joule  heat).  As  a  result,  the  temperature  of 
the  carbon  material  increases,  causing  desorption. 

The  adsorption  process  with  the  reversed  gas  flow 
(inlet  through  the  outer  annular  pipe  and  outlet  through  the 
central  pipe),  which  is  also  feasible  in  this  geometry,  is  also 
considered. 

Two  different  models  of  the  adsorber  will  be  present.  In  the  first  one,  the  adsorbent 
bed  is  regarded  as  homogeneous.  In  the  second  one,  a  more  realistic  picture  of  the  adsorbent 
bed  is  used,  taking  into  account  its  structure,  made  of  several  layers  of  activated  carbon  cloth. 
In  both  models,  uniform  conditions  within  the  adsorbent  particles  (ACFC  fibers)  are  assumed, 
i.e.  both  models  can  be  regarded  as  macroscopic  [4], 

Both  mathematical  models  are  obtained  as  sets  of  coupled  nonlinear  partial  and 
ordinary  differential  equations,  defining  the  mass  and  heat  balances  in  the  system.  The  models 
also  include  expressions  for  calculation  of  the  equilibrium  and  transport  parameters. 

The  model  equations  are  solved  numerically.  The  PDEs  are  first  approximated  by 
larger  sets  of  ordinary  differential  equations,  using  the  method  of  orthogonal  collocations  [5], 
The  resulting  sets  of  ODEs  are  further  solved  numerically,  leading  to  the  concentration  and 
temperature  profiles  in  time  and  space  in  the  adsorber. 

The  obtained  models  are  used  for  simulation  of  the  processes  of  adsorption, 
electroresistive  heating  and  electrothermal  desorption,  as  well  as  consecutive  adsorption- 
desorption,  for  a  particular  system. 


f  GAS  IN 


Figure  3.  A  schematic  representation  of  the 
annular  -  radial  flow  -  ACFC  adsorbent  bed 
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3.  Mathematical  model  of  the  single-cartridge  adsorber  with  homogeneous  adsorbent  bed 

In  postulating  this  model  it  was  assumed  that  the  layers  of  the  ACFC  were  coiled  so  tightly, 
that  the  contact  between  them  was  very  intimate  and  the  adsorbent  density  was  uniform 
throughout  the  bed,  so  that  the  whole  adsorbent  bed  could  be  treated  as  homogeneous. 

The  following  additional  assumptions  were  used  in  setting  up  this  model: 

•  The  mass  and  heat  transfer  resistances  on  the  particle  scale  can  be  neglected. 

•  There  are  no  changes  of  the  concentrations  and  temperatures  in  the  axial  direction,  i.e. 
only  the  gradients  in  the  radial  direction  have  to  be  taken  into  account,  so  that  one¬ 
dimensional  model  can  be  used  [4], 

•  Initially,  the  adsorbate  concentrations  in  both  phases  are  in  equilibrium  and  uniform 
throughout  the  adsorbent  bed  (for  the  cases  of  adsorption  and  electroresistive  heating 
without  desorption,  they  are  both  0).  The  temperatures  of  both  phases  are  initially 
equal  and  uniform  throughout  the  adsorbent  bed. 

•  The  fluid  phase  is  an  ideal  gas  mixture  of  the  inert  and  the  adsorbate  (for  the  case  of 
electroresistive  heating  without  desorption,  it  is  pure  inert). 

•  The  gas  pressure  is  constant  and  uniform  throughout  the  bed. 

•  Heat  dispersion  in  the  radial  direction  of  the  adsorbent  bed,  both  in  the  solid  and  in  the 
fluid  phase,  is  taken  into  account. 

•  The  gas  in  the  inlet  tube  and  in  the  space  around  the  adsorbent  bed  is  ideally  mixed. 

•  All  physical  parameters  and  coefficients  are  considered  as  constants. 

•  The  electric  resistivity  of  the  ACFC  adsorbens  is  temperature  dependent. 

•  The  electric  power  is  supplied  under  the  constant  voltage  conditions  during 
electroresistive  heating  and  electrothermal  desorption. 


3.1.  Mass  and  energy  balances 

Based  on  these  assumptions,  the  mathematical 
model  of  the  annular  radial-flow  adsorption  bed 
presented  in  Fig.  3  was  obtained  by  writing  the  material 
and  energy  balances  for  a  differential  element  of  the 
cartridge  volume,  as  shown  in  Figure  4.  By  letting  the 
thickness  of  this  element  A r  to  zero,  these  balances  are 
obtained  as  the  following  partial  differential  equations: 


Adsorbate  balance  for  the  solid  phase  within  the 
adsorbent  bed: 

p.lH^XC-c')  a) 

Ot 


Figure  4  A  differential  volume  element,  used 
for  modeling 
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Adsorbate  balance  for  the  gas  phase  within  the  adsorbent  bed: 


d£  G  dC 
dt  (2rnHeb)  dr 


+  (kma)(C-C*)  = 


D  d 


r  dr 


dC_ 

dr 


Heat  balance  for  the  solid  phase  within  the  adsorbent  bed: 


Pi 


dt  L 


\cp,+cplc/)T,]  =  ^+AHmhpt^-+ 

J  dV  dt  r  or 


dq  D'f  d 


dT 


r — - 

K  dr  j 


-( ha)(Ts-T ) 


(2) 


(3) 


Heat  balance  for  the  gas  phase  within  the  adsorbent  bed: 

<Kc„+cwQ  dTs 

- I  if.  -i-  i  ii#  i  — - 

gdtl 


d  f/ 

P*— L<' 


c  +  c 

Pg  pv 


c)r„]  =  -- 


2nzHzh 


dr 


GC,Js  5C  +  (ha)(T-T)+D’S  d 


2rnHsb  dr 


r  dr 


r  dT  \ 
r — - 
v  dr  J 


(4) 


The  boundary  conditions  for  equations  (1-4)  are: 


r=r\\  Z) 


dC_ 
‘  dr 


=  c-cit,  d; 


hs 


r=rx 


st. 

dr 


=  ksl(T,-T„),  D* 


dT 


dr 


=  K^g-TgU) 


n  dC\ 
r=r2:  ~Dm—\ 

or 


hs  df 


=  C-C  -D  — — 

^  ^ot  ’  ~ 

or 


=  -d, 


M  dTt 


dr 


(5) 


=  kJT!~T„)  (6) 


Using  the  assumption  of  ideal  mixing  in  the  central  inlet  tube  and  in  the  annular  space 
around  the  adsorption  bed  (the  outlet  tube),  the  mass  and  heat  balances  for  the  gas  phase  in 
these  two  regions  are  obtained  in  the  form  of  ordinary  differential  equation: 

Adsorbate  balance  for  the  gas  in  the  inlet  tube: 


dC, 


nH9g  —±-  =  GCin  -  GCit  +  kml  (2rlnH)(l  -  sb )(C*  _  -  C„ )  +  2rxiiH&b (C^  -  C, ) 


(7) 


(8) 


dt 

Heat  balance  for  the  gas  in  the  inlet  tube: 

jX(Cp‘  +C'"C">r*'']  =  G(C«  +C,»C»>7ir.  "G<C«  +C„.CI,K , 

+K  )(1  -  ■=»  )(?;H  -  Tr, )  +  hsl  (2 r^Hs,  )(rsH  -  7V, ) 

Adsorbate  balance  for  the  gas  in  the  outlet  annular  tube: 

(ri  ~r2  )%H Pg  —fr  =  GC\  -GC„+*.!(2r2nffXl-etXC,|  -CJ  +  2rinHh(C\r_r  -CJ 

Heat  balance  for  the  gas  in  the  outlet  annular  tube: 

(r?  -r2)nHp  —\(c  +c  C  AT  ,~\  =  G(c  +c  C\  )T ,  -G(c  +c  CAT 

V3  2/  r  g  ,  |_V  pg  pv  old  got  J  V  pg  pv  |  r=r2J  g\r=r2  v  Pg  pv  ot  A  g 


L v  Pg  Pv  ol '  g°l  J  v  Pg  Pv  I  r=r2'  g\r=r2  v  pg  pv  vis  got  (10) 

+hs2(2r2nH)(l-eb)(T  -Tgot)  + hg2(2r2nHsb)(T,  -Tgol)-hwg(2r2nH)(Tgot-Ta) 


(9) 
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(11) 


The  initial  conditions  for  equations  (1-4)  and  (7-10)  are: 

tO,  re[rur2]:  Tg=Ts=Tgit=Tgot=Tgin=Tgo=Ta=Tp,  C=Cit=Cot=Cul=C*=(\Hq)=Cp 


The  notations  of  the  main  variables  used  in  Figure  4  and  equations  (1-11)  are:  C  - 
adsorbate  concentration  in  the  gas  phase,  q  -  adsorbate  concentration  in  the  solid  phase,  C*  - 
adsorbate  concentration  in  the  gas  phase  in  equilibrium  with  the  solid  phase,  Tg  -  gas 
temperature,  Ts  -  solid  temperature,  t  -  time,  r  -  the  radial  position  in  the  cartridge,  r\  and  r2  - 
the  inner  and  the  outer  diameter  of  the  radial  adsorbent  bed,  r 3  the  diameter  of  the  adsorber,  H 
-  the  cartridge  height,  8b  -  bed  porosity,  G  -  molar  flow-rate  of  the  gas  inert.  The  subscript  in 
denotes  the  inlet,  o  the  outlet,  it  the  inlet  tube,  ot  the  outlet  annular  tube,  a  the  ambient  and  p 
the  previous  (initial)  state.  6Qej  is  the  heat  Joule  heat  generated  in  the  differential  element  of 
volume.  A  complete  list  of  the  notations,  including  the  physical  parameters  and  coefficients 
used  in  these  equations,  can  be  found  in  the  nomenclature,  at  the  end  of  this  report. 

These  model  equations  are  valid  for  all  three  processes:  adsorption,  electroresistive 
heating  and  electrothermal  desorption,  assuming  the  proper  definition  of  the  variables  (e.g.  for 
adsorption  8Qei= 0,  Cp= 0  and  C,„=conste0,  for  electroresistive  heating  Cp= 0  and  Cin= 0,  while 
for  electrothermal  desorption  Cp=const^0  and  C,„=0). 


3.2.  Additional  equations 

For  a  complete  mathematical  model  besides  the  balance  equations  and  their  initial 
and  boundary  conditions,  it  is  necessary  to  define  relations  for  calculating  the  equilibrium  and 
transport  parameters.  A  selection  of  these  relations,  based  on  literature,  is  given  below. 


Calculation  of  the  equilibrium  parameters  (adsorption  isotherm) 

The  equilibrium  relation  C*=0(q )  is  one  of  the  basic  data  needed  for  correct  modeling 
of  adsorption  systems.  In  the  presented  model  the  Dubinin-Radushkevich  equation  was  used, 
as  it  had  been  shown  to  describe  reliably  adsorption  of  organic  vapors  on  ACFC  [3]: 


f 

W  =  W()  exp 

v 


R.T 


Tn^ 
Pa 


W 


JJ 


(12) 


In  equation  (13)  IF  is  the  volume  of  the  adsorbed  phase  per  unite  mass  of  adsorbent, 
Wo-  the  volume  of  the  micropores  per  unit  mass  of  adsorbent,  E-  the  adsorbate  energy  of 
adsorption,  />°-the  adsorbate  saturation  vapor  pressure,  p,\  -  the  adsorbate  partial  pressure,  T— 
temperature  and  Rg  -  the  gas  constant.  The  vapor  pressure  was  calculated  using  the  Wagner 
equation: 


In 


(  0  \ 

p  ' 

'  VPAx  +  VPBx15  +  VPcx3  +  VPDx6 '' 

UcJ 

l  J 

(13) 


where  x  =  1  -  (Ts  /Tc) ,  pc  and  Tc  are  the  critical  pressure  and  temperature  and  VP  a,  VPb,  VPc 
and  VPd  are  Wagner  constants. 
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Based  on  equation  (12),  the  following  relation  between  the  equilibrium  adsorbate 
concentration  in  the  gas  phase  C*  and  its  concentration  in  the  solid  phase  q  and  the 
temperature  Ts,  was  obtained: 


C  = 


p  exp 


In 


P  aK. 


(14) 


(Mj  and  p a  are  the  adsorbate  molar  mass  and  density). 

Depending  on  the  particular  system,  another  adsorption  isotherm  relation  can  be  used,  either 
based  on  theory  or  on  experimental  results. 


Calculation  of  the  transport  parameters 

The  calculation  of  the  main  transport  parameters  was  based  on  different  correlations 
published  in  the  literature.  Most  of  the  correlations  used  in  this  work  are  general,  and  do  not 
correspond  particularly  to  the  fibrous  structure  of  the  adsorption  bed.  They  could  be  replaced 
with  more  specific  correlations,  if  they  become  available.  Here  are  the  procedures  for 
calculation  of  the  main  transport  parameters  used  in  the  model: 


•  Heat  transfer  coefficient  between  the  ACFC  and  the  gas  in  the  inlet  or  outlet  tube  hsj 
and  hS2 

The  criterial  equation  suggested  in  Ref.  [13],  was  used  for  calculation  of  this  heat 
transfer  coefficient: 


_  NuHk 

ns  1/2  — 


Nuh  =i 


H 


1/6 


0.825+  °-387  RaH 


RaH  — 


[l  +  (0.492  /  Pr)9/16  ] 


Pr  = 


a  = 


av 

P(Cre+^vC) 

k 

k 

p  Acpg  +CFC) 


(15) 

(16) 

(17) 

(18) 
(19) 


In  these  equations,  Nuh  is  the  Nuselt  numper,  Pr  -  the  Prantl  number,  Ran  -  the 
Rayleigh  number,  k  -  the  thermal  conductivity  of  the  gas  phase,  g  -  the  gravity  constant,  P- 
the  reciprocal  value  of  the  temperature  of  the  film  (approximately  calculated  as  an  arithmetic 
mean  of  the  temperatures  of  the  solid  surface  Ts  and  the  bulk  temperature  of  the  gas  77),  a  - 
the  thermal  diffusivity  and  v  and  p  -  the  kinematic  and  the  dynamic  viscosity  of  the  gas  phase, 
respectively. 
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•  Radial  thermal  conductivity  of  the  gas  phase  within  the  adsorbent  bed  Dthg 

A  correlation  proposed  by  Edwards  and  Richardson  [7]  was  used: 

1  _  Dgf  _  0.73Sfr  0.5 

Pe  PgvDp  (cpg  +  cpvC)  RePr  \  +  (9.1zb)l{RePr) 


(20) 


Re  = 


E 

G 

2rnHpg 


(21) 

(22) 


where  Pe  is  the  Peclet  number,  Re  -  the  Reynolds  number,  Dp  -  the  particle  diameter  and  v  - 
the  superficial  gas  velocity.  The  value  of  the  Pe  number  changes  with  time  and  the  radial 
position  in  the  bed,  while  the  superficial  velocity,  and  consecutively,  the  Re  number  changes 
with  the  radial  position  in  the  bed.  Accordingly,  the  radial  thermal  diffusivity  is  a  function  of 
both  time  and  the  radial  position  in  the  adsorbent  bed. 


•  Gas  to  solid  mass  transfer  coefficient  in  the  ACFC  adsorbent  bed  km 
This  coefficient  was  calculated  using  the  following  criterial  equation  [6]: 


Sh  =  0.9  IT  Re0'49  Sc1/3 

(23) 

Sc  =  » 

Pg^AB 

(24) 

k  D 

5  1 

II 

$5 

(25) 

d<b 

where  Sh  is  the  Sherwood  number,  Sc  -  the  Schmidt  number,  Dab  -  the  molecular  diffusivity 
of  the  adsorbate  in  the  gas  inert  and  ¥  -  the  shape  factor  (0.81  for  cylindrical  particles). 

The  molecular  diffusivity  DAb  was  calculated  using  the  Chapman-Enskog  correlation 

[6]: 


0.001858T372  | 

ma  +  mb 

*  i 

2MaMb 

P^AB&D 


(26) 


with:  a  AB  =  —  +  —  -  the  characteristic  length  of  the  binary  pair 

ct  =2.44  ( TJPc)m  -  the  characteristic  length  of  a  single  molecule 
Q,  -(44.54J~~  +  1.91  ir-- )“.  ,he  comsion  inte^al 

T'  =  kbTg/tJB 


8 AB  ~  (8A8b) 

Sa  -  the  characteristic  Lennard  -Jones  energy 
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s/k  =  0J5Tc 

ho  -  Boltzmann’s  constant 

•  Gas  to  solid  heat  transfer  coefficient  in  the  ACFC  adsorbent  bed  h: 

Based  on  the  analogy  between  mass  and  heat  transfer  in  the  adsorbent  bed  [6,  9],  a 
criterial  equation  analogous  to  equation  (25)  was  used: 

Nu  =  0.9  IT  Re0  49  Pr1/3  (27) 


•  Radial  mass  dispersion  coefficient  in  the  gas  phase  Dm 

Based  on  the  analogy  between  mass  and  heat  transfer  [9]  and  the  equation  for  the  heat 
dispersion  coefficient,  the  following  was  used: 


D„  = 


Dhg 


(c  +c  C) 

V  pg  pv^S 


(28) 


Calculation  of  the  Joule  heat 

An  important  term,  differentiating  electroresistive  heating  from  other  processes,  is 
'bQeildV,  in  equation  (3).  It  represents  the  Joule  heat  generated  in  a  unit  bed  volume.  For  the 
case  of  constant  electric  voltage  supply,  this  term  can  be  derived  through  the  following  steps: 

•  Calculation  of  the  specific  electric  resistivity  of  ACFC,  which  is  temperature 
dependent.  A  linear  relationship,  suggested  in  Ref.  [3]  was  used: 

p  =  f(Ts )  =  (p0  (1  +  b(Ts  -  Tr  )))  =  f(r,  t )  (29) 

where  po  is  the  electric  resistivity  at  the  referent  temperature  Tr  and  b  is  the  temperature 
coefficient  of  the  electric  resistivity. 

•  Calculation  of  the  electric  resistance  of  the  differential  volume  element  defined  in 
¥ig.4(dV  =  2rnHdr): 


dRel(r,t) 


P(T,)H 

dA 


P(T,)H 

2rndr 


(30) 


•  Calculation  of  the  electric  current  through  the  differential  volume  element  dV,  under 
constant  voltage  conditions: 


dl  (r,  t ) 


U 

dRel(r,t) 


2  rnUdr 
P  (Ts)H 


(31) 


•  Calculation  of  the  Joule  heat  generated  in  the  differential  volume  element  dV : 


80- 


Udl(r,t ) 


2nzU2dr 
P  (T,)H 


(32) 


•  Calculation  of  the  heat  generated  per  unite  volume: 

80„  U 2 

dV  p(Ts)H2 


(33) 
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The  electric  resistivity  p(Tv)  is  calculated  for  the  local  value  of  the  solid  phase 
temperature.  As  this  temperature  varies  both  with  time  and  with  the  radial  position  in  the  bed, 
so  does  the  Joule  heat  supplied  to  the  element  of  the  volume  in  question. 


3.3.  Numerical  solution  of  the  model  equations 

The  first  step  of  the  numerical  solution  was 
approximation  of  the  partial  differential  equations  (1-4) 
by  a  larger  set  of  ordinary  differential  equations,  by 
applying  the  method  of  orthogonal  collocation  [5].  Before 
applying  this  method,  the  radial  coordinate  of  the 
adsorption  bed  was  transformed  into  a  nondimensional 
one,  using  the  following  relation: 


u  = 


(34) 


Figure  5  gives  a  schematic  representation  of  the 
distribution  of  the  collocation  points  in  the  bed.  N  internal 
collocation  points  were  used,  plus  the  border  points  u= 0 
(r=r\)  and  u=  1  (r=r?). 

Here  are  the  resulting  equations: 


Figure  5  Schematic  representation  of  the 
distribution  of  the  collocation  points 


^  =  ^(c,-c;),  ;=o,...^v+i 

dt  pb 


N+ 1 


dC\ 


c, 


j= 0 


kci 


(35) 

D  V  B  C 

m  L-t  i,J  J 


N+ 1 


dt  27i//s,o  A(Am,  +/,)  p 


(Ct-Ct)  +  - 


i= o 


(Am  +  rx)p 


i=l,...,N+l 


(36) 


dT„ 


6  a 


dt  dVpb(cps  +  cplq^) 


D, 


hs 


N+ 1 


P  biCps+Cpldi) 


(uA  +  rt)ph(ci)S+cplq)  /=0 


ZBi,lTv’  (37> 


i=l . N 


N+ 1 


dTg  (to)(r„-r8,) 


7=0 


7V+1 


Gc  T  Y  AC 

pv  gi  Lu  IJ^J 
7=0 


*  Pg(Sg  +SvQ)  27tffe*p,,A(Kl.A  +  r1)  27cff6^A(Ml.A  +  r1)(cwr  +  c^Cf) 


pv 

T  dCt 

c  T  .  — L 

dt 


(38) 


-  +  ■ 


7V+1 


(C*  +  CPvCi )  (U,A  +  h  )P*  (Cpg  +  CpvCi ) 


Ev«.  i=1-iV 


7=0 


The  mass  and  heat  balances  for  the  inlet  and  outlet  gas  tubes  were  also  rearranged: 

7C,  G(C„  -C,)  {  &ml (2^71/7 )(1  - )(Cq  - Cit )  {  2r17Iffe,(C0  -C,) 
dt  rfnHpg  rxnHpg  rl2nHpg 


(39) 
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Jrg,  G(c„+cpfCJTgb,  GT ’?1  |  hsl  (2r17i//)(l  -  eb  )(7^0  -  7A, )  < 

dt  r;KHpg(cpg+cpvC,l)  tfnH  pg  /y7i//pg(cyg +c/n,C,() 

\i  (2r17t//sf,  )(Tg0  -  7A, )  CpJju  dCjt 

K^H pg  (cpg  +  cpvCit )  (c/)g  +  cpvCit )  dt 

dC o,  GjC^  -CJ  {  kj,  2r27r//)(l-s,)(C;+1  -CJ  {  2r27i//s,(C;v+1  -CJ 

dt  (tf  - r; )nH pg  (r32-r22)ntfpg  (r32-r22)^pg 


dT 


got  _ 


G(Cpg  +  C  pVC N+l^TgN+l 


GT. 


got 


dt 


,  42(2r27i//)(l-8j(^+1-7;j  | 


(r3  -  r2  )7i/f pg  (c  +  c  Cc< )  (r3  -  r2  )7tffpg  (r3  -  r2  )7ttf pg  (c  +  cCot ) 


hgl(2r2KHzb)(TgN+x  -T  )  h(2r,nH){Tgot -Ta) 


rri 

^ pv  got 


dC„. 


(42) 


(T  ~  G  )nH pg  (cpg  +  cpvCot )  (r3  -  r2  )7iffpg  (c/)g  +  cpvCot )  (cpg  +  cpvCot )  dt 

The  boundary  conditions  (6-7)  are  also  transformed  into  the  following  algebraic 
expressions: 

N+\  i  N+ 1 


Q  = 


1  TV+1 

A  +  A.tI-VA 

_ A  7=1 _ 

i  • 


T  = 


-j  TV+1  i  TV+1 

+ a**  x  1 4  a  *„t;« + A"  -r  z  Aoj,i 

— .  a,= - ^ — (43> 


A.+A^A. 


a,2+A*'^A,. 


A-A.tZ-W’.,- 


c 

^N+ 1 


A  7=0 


T  = 

?  ^  sTV+1 


1  +  Z> 


1  TV  1  TV 

-  A*’  7IW,  -A%  tEU 

_  T  = 

5  ^gTV+1 


W+l, TV+1 


/,  +  J_  // 

"s2  T  ^  ''W+l.JV+l 


h,,-,  +  D,/,g  -  7 


5  2 


t  ^  "*TV+1,TV+1 


Aij  and  T?y  are  the  first  and  second  order  derivatives  of  the  Lagrangian  polynomials  /, 
used  in  the  method  of  orthogonal  collocations  [5]: 


Aj  = 


f  dlj  } 

.  du  , 

V  /  u=Uj 


u,A  +  rx 
Bu=-^- 


v  Jl/2  , 

V  s  u=ui 


+  - 


1 


r  dl  A 


v  j 

v  y  u=u: 


(45) 


(A  =r2-ri  is  the  bed  thickness.) 

Equations  (35-42)  represent  a  set  of  4 AM-6  nonlinear  ODEs  which,  together  with  6 
algebraic  equations,  define  the  concentrations  and  temperatures  in  the  solid  and  in  the  gas 
phase  in  all  collocation  points  of  the  annular  adsorption  bed,  as  well  as  in  the  gas  phase  in  the 
inlet  and  in  the  outlet  tube.  This  set  of  equations  was  solved  numerically,  using  the  Matlab 
function  odel5s. 

The  derived  model  was  used  for  simulation  of  the  concentration  and  temperature 
profiles  of  both  phases  in  different  points  of  the  ACFC,  as  well  as  in  the  inlet  and  outlet  tube, 
for  3  different  processes:  adsorption  on  an  initially  clean  ACFC  bed,  electroresistive  heating 
of  a  clean  ACFC  bed  and  electrothermal  desorption  from  a  previously  saturated  ACFC  bed. 
The  parameters  used  for  simulation  are  listed  in  Table  1.  They  correspond  to  the  experimental 
system  defined  in  Ref.  [3]:  methyl  ethyl  ketone  (MEK)  as  adsorbate,  American  Kynol  ACC- 
5092-20  as  adsorbent  and  nitrogen  as  carrier  gas.  The  equilibrium  and  most  of  the  transport 
parameters  were  calculated  using  the  expressions  given  in  Section  3.2,  while  for  the  rest, 
arbitrary  values,  given  in  Table  1,  were  used. 


(44) 
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Table  1.  Parameter  values  used  for  simulation  of  adsorption,  electroresistive  heating  and  electrothermal 
desorption,  using  the  model  of  a  homogeneous  adsorption  bed 


Variable 

Notation  and  value 

Inner  diameter  of  the  cartridge 

r^O.95  cm 

Outer  diameter  of  the  outlet  annular  tube 

r3=3.55  cm 

Bed  thickness 

A=r2-ri=0.55  cm 

Bed  axial  dimension 

H=  30  cm 

Bed  porosity 

8/, =0.72 

Particle  (fiber)  diameter 

Z)p=13xl0"4  cm 

Specific  surface  area 

a=13.65  cm2/cm3 

Gas  pressure 

£>=101325  Pa 

Inlet  adsorbate  concentration  in  the  gas 
phase  during  adsorption 

Cin= 0.001  mol/mol 

Inlet  adsorbate  concentration  in  the  gas  phase  during  desorption 

Cin= 0  mol/mol 

Initial  adsorbate  concentration  in  the  solid  phase  for  desorption 

gp=0.003  mol/g 

Inlet  gas  temperature 

rem=293.15  K 

Ambient  temperature 

Ta=293.15  K 

Mass  transfer  coefficient  in  bed 

A:, ,,=0.0003  mol/(cm2s) 

Mass  transfer  coefficient  from  the  solid  phase  to  the  gas  phase  in 
the  inlet  tube 

kml=0.0003  mol/(cm2s) 

Mass  transfer  coefficient  from  the  solid  phase  to  the  gas  phase  in 
the  outlet  annular  tube 

km2= 0.0003  mol/(cm2s) 

Solid  to  gas  heat  transfer  coefficient 

h= 9x1  O'4  W/(cm2K) 

Heat  transfer  coefficient  from  the  gas  phase  in  bed  to  the  gas 
phase  in  the  inlet  tube 

hgl- 8x1 0-3  W/(cm2K 

Heat  transfer  coefficient  from  the  gas  phase  in  bed  to  the  gas 
phase  in  the  outlet  annular  tube 

hg2= 2x1 0'3  W/(cm2K 

Gas  to  ambient  heat  transfer  coefficient  (heat  losses) 

hwe= 5x1 0-6  W/(cm2K) 

Heat  diffusivity  of  the  solid  phase 

Z)/,s=  1.2x1 0‘4  W/(cmK) 

Heat  of  adsorption 

AHlld, =66200  J/mol 

Adsorbent  bed  density 

p*=0.221  g/cm3 

Density  of  the  inert  gas 

p„=0.000416  mol/cm3 

Adsorbate  density 

p,4=0.81  g/cm3 

Adsorbate  molar  mass 

Ma=12.\01  g/mol 

Inert  gas  molar  mass  (nitrogen) 

Mg=28.02  g/mol 

Specific  heat  capacity  of  liquid  adsorbate 

c0/=157.9  J/(molK) 

Specific  heat  capacity  of  the  inert  gas 

cve= 29.13  J/(molK) 

Heat  capacity  of  the  solid  phase 

cps=0.7l  J/(gK) 

Specific  heat  capacity  of  gaseous  adsorbate 

CpV=  100.9  J/(molK) 

Dynamic  viscosity  of  the  inert  gas 

p=1.78xl0'4  g/cm/s 

Inert  gas  thermal  conductivity 

&=2.6xl0'4  W/cm/K 

Referent  temperature 

Tfi=293.15  K 

Critical  temperature  of  the  adsorbate 

TcA= 536.8  K 

Critical  temperature  of  the  inert  gas  (nitrogen) 

TcB=  126.2  K 

Critical  pressure  of  the  adsorbate 

pcA=4.26xlO+6  Pa 

Critical  pressure  of  the  inert  gas  (nitrogen) 

£>cS=3.39x10+6  Pa 
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Wagner  constants  for  the  adsorbate 

VPA=-1 .1X41 6 

VPB=\.l\i)6\ 

FPc=-3.6877 

FPfl=-0.75169 

Total  volume  of  micropores  (D-R  equation) 

1F0=0.748  cm3/g 

Gas  constant 

i?g=8.314  J/(molK) 

Gravitation  constant 

g-=9.81xl02  cm/s2 

Boltzmann’s  constant 

4=1.38048xl0'23  J/K 

Adsorption  energy  of  the  adsorbate  (D-R  eq.) 

£=14.43xl0+3  J/mol 

Electric  resistivity  at  referent  temperature  TR 

p0=0.202  Qcm 

Temperature  coefficient  of  the  bed  electrical  resistivity 

6=1.02e'5  1/T 

3.4.  Simulation  of  adsorption 

An  important  step  of  the  TSA  process  is  the  adsorption  part.  Using  the  model 
presented  here,  simulation  of  adsorption  on  a  previously  clean  ACFC  bed,  for  constant  inlet 
gas  concentration  and  temperature,  was  performed.  For  this  case,  the  model  equations  were 
solved  for  the  following  conditions:  S£>e/=0,  Cp= 0,  rp=const,  C,„=const,  7’;„=const. 

Some  of  the  simulation  results  obtained  for  the  parameter  values  given  in  Table  1 
and  for  10  collocation  points  (A+2=10)  are  shown  in  Figures  6-9. 

Figures  6  and  7  show  the  space  and  time  distributions  of  the  concentrations  (Fig.  6) 
and  temperatures  (Fig.  7)  in  the  solid  (figs  (a))  and  in  the  gas  phase  (figs  (b)).  In  these  figures, 
lines  of  different  colors  correspond  to  different  collocation  points  in  the  bed,  the  distribution 
of  which  is  shown  in  the  figures  as  well,  by  colored  dots  (e.g.  the  red  lines  and  red  dots 
correspond  to  the  inner  boundary  point  r=r\,  while  the  gray  ones  correspond  to  the  outer 
boundary  point  r=r2).  The  concentrations  and  temperatures  of  the  gas  phase  in  the  inlet  and  in 
the  outlet  tube  are  also  shown,  as  dashed  black  lines.  For  the  parameters  used,  these  lines 
overlap  with  the  border  lines  corresponding  to  the  gas  phase  in  the  adsorption  bed. 

Figures  6  and  7  actually  show  the  propagation  of  the  concentration  and  temperature 
fronts  in  the  radial  direction  of  the  bed. 


Figure  6  Simulation  of  adsorption  -  time  profiles  of  the  concentrations  in  10  collocation  points  in  the  ads.  bed  (colored  lines) 
and  in  the  inlet  and  outlet  tubes  ( — ),  for  G=0.1  mol/s  :  a)  Solid  phase;  b)  Gas  phase 
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Figure  7  Simulation  of  adsorption  -  time  profiles  of  the  tempaeratures  in  10  collocation  points  in  the  ads.  bed  (colored  lines) 
and  in  the  inlet  and  outlet  tubes  ( — ),  for  G= 0. 1  mol/s  :  a)  Solid  phase;  b)  Gas  phase 


The  simulation  was  also  used  to  investigate  the  influence  of  the  inert  gas  flow-rate  on 
the  adsorption  results.  The  simulated  concentrations  in  the  solid  phase  at  the  outer  surface  of 
the  adsorption  bed  (for  r=r2)  and  in  the  outlet  gas  streams,  obtained  for  3  different  flow-rates, 
are  shown  in  Figures  8(a)  and  8(b),  and  the  corresponding  temperatures,  in  Figures  9(a)  and 
9(b). 
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Figure  8.  Simulation  of  adsorption  -  influence  of  the  gas  flow-rate  on  the  concentrations: 
a)  Concentration  in  the  solid  phase  at  r=r2 ,  b)  Outlet  gas  concentration 


It  has  to  be  mentioned  that  the  expressions  for  calculation  of  the  mass  and  heat 
transfer  parameters  km  and  h,  presented  in  Section  3.2.  could  not  be  used  in  these  simulations, 
as  they  gave  values  which  caused  serious  convergence  problems.  The  values  given  in  Table  1 
were  used  instead.  The  calculated  values  of  h  were  used  only  for  simulation  of  electroresistive 
heating. 
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Figure  9.  Simulation  of  adsorption  -  influence  of  the  gas  flow-rate  on  the  temperatures: 
a)  Temperature  of  the  solid  at  r=r2 ,  b)  Outlet  gas  temperature 


3.5.  Simulation  of  electroresistive  heating 

In  order  to  get  a  reliable  description  of  the  electrothermal  desorption  process,  it  is 
important  to  correctly  model  pure  electroresistive  heating  of  an  adsorbate-free  bed,  first.  This 
model  would  enable  to  predict  the  temperature  profiles  in  the  bed.  Having  some  experimental 
measurements  of  the  temperatures  would  enable  estimation  of  the  parameters  defining  the 
heat  generation  and  heat  transfer  separately,  and  using  them  in  modeling  of  electrothermal 
desorption. 

The  model  presented  in  Section  3.1  can  as  well  be  used  for  electroresistive  heating, 
assuming  that  Joule  heat  is  generated  inside  the  solid  phase  and  that  both  the  gas  and  the  solid 
phase  are  adsorbate-free  at  all  times:  0,  Cp= 0,  7),=const,  Q„= 0,  7/„=const.  Accordingly, 

all  material  balance  equations  and  their  appropriate  boundary  conditions  can  be  omitted,  and 
the  solid  and  gas  heat  balances  for  the  adsorption  bed,  together  with  the  heat  balances  for  the 
inlet  and  the  outlet  tube,  can  be  solved  separately.  After  applying  the  method  of  orthogonal 
collocations,  the  model  reduces  to  2N+2  ODEs  and  4  algebraic  equations. 

The  simulation  results  presented  here  (in  Figures  10-12)  were  obtained  for  constant 
electric  voltage  supply.  Ten  collocation  points  (7V+2=10)  were  used  for  numerical  solution  of 
the  PDEs.  The  simulated  time  profiles  of  the  solid  and  gas  temperatures  in  different 
collocation  points  in  the  bed  are  shown  in  Figure  10,  as  colored  lines,  and  the  gas 
temperatures  in  the  inlet  and  in  the  outlet  tube  as  black  dashed  lines.  The  influence  of  the  gas 
flow-rate  is  shown  in  Figure  11,  and  the  influence  of  the  applied  electrical  voltage  in  Figure 
12.  In  these  two  figures,  the  simulated  time  profiles  of  the  outlet  gas  temperature  and  the  solid 
temperature  at  the  outer  surface  of  the  cartridge  are  shown. 
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(a)  (b) 


Figure  10  Simulation  of  electroresistive  heating  -  time  profiles  of  the  temperatures  in  10  collocation  points  in  the  ads.  bed 
(colored  lines)  and  in  the  inlet  and  outlet  tubes  ( — ),  for  G=0.02  mol/s  and  U=  12  V  :  a)  solid  phase,  b)  gas  phase 


(a)  (b) 


t  (s)  t  (s) 

Figure  1 1  Simulation  of  electroresistive  heating  -  influence  of  the  gas  flow-rate  on  the  temperatures  for  U=12  V: 
a)  Temperature  of  the  solid  at  r=r2,  b)  Outlet  gas  temperature 
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Figure 


12  Simulation  of  electroresistive  heating  -  influence  of  the  electric  voltage  on  the  temperatures  for  G=0.02  mol/s: 
a)  Temperature  of  the  solid  at  r=r2,  b)  Outlet  gas  temperature 
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3.6.  Simulation  of  electrothermal  desorption 

In  electrothermal  desorption,  a  previously  adsorbed  adsorbate  is  desorbed  owing  to 
Joule  heat  generated  inside  the  adsorbent  material  and  a  flow  of  pure  inert  through  the 
adsorbent  bed.  The  model  developed  in  Section  3.1  was  used  for  simulation  of  electrothermal 
desorption,  for  the  case  of  thermal  and  concentration  equilibrium  and  uniform  concentration 
and  temperature  distributions  in  the  adsorbent  bed  in  the  initial  state,  and  for  constant  electric 
voltage  supply  during  electrothermal  desorption  ( Cp=consteO, Theorist,  Cin= 0,  T,,„=const  and 
f/=const). 

Some  of  the  simulation  results  obtained  with  10  collocation  points  are  shown  in  13-16. 
The  time  profiles  in  different  collocation  points  in  the  adsorption  bed  and  in  the  inlet  and 
outlet  tubes)  are  shown  in  Figure  13  (concentrations  q  and  Q  and  Figure  14  (temperatures  Ts 
and  Tg ).  Figure  15  illustrates  the  influence  of  the  gas  flow-rate,  and  Figure  16  the  influence  of 
the  applied  voltage  on  the  concentrations  and  temperatures  corresponding  to  the  outer 
boundary  of  the  adsorbent  bed  and  the  outlet  gas  stream. 
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Figure  13  Simulation  of  electrothermal  desorption  -  time  profiles  of  the  concentrations  in  10  collocation  points  in  the  ads. 
bed  (colored  lines)  and  in  the  inlet  and  outlet  tubes  ( — ),  for  G=0.03  mol/s  and  U=  10  V:  a)  solid  phase;  b)  gas  phase 
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Figure  14  Simulation  of  electrothermal  desorption  -  time  profiles  of  the  temperature  in  10  collocation  points  in  the  ads.  bed 
(colored  lines)  and  in  the  inlet  and  outlet  tubes  ( — ),  for  G=0.03  mol/s  and  G=10  V:  a)  solid  phase;  b)  gas  phase 
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t  (s)  t  (s) 

Figure  15  Simulation  of  electrothermal  desorption  -  influence  of  the  electric  voltage  on  concentrations  and  temperatures  at 
r=r2  (solid)  and  in  the  outlet  tube  (gas)  for  G=0.03  mol/s:  a)  Concentrations  and  b)  Temperatures 


t  (s)  t  (s) 

Figure  16  Simulation  of  electrothermal  desorption  -  influence  of  the  electric  voltage  on  concentrations  and  temperatures  at 
r=r2  (solid)  and  in  the  outlet  tube  (gas)  for  G=0.03  mol/s:  a)  Concentrations  and  b)  Temperatures 


3.7  Simulation  of  consecutive  adsorption-desorption 

The  TSA  process  with  electrothermal  desorption  step  consists  of  consecutive 
adsorption  and  desorption  steps.  The  model  defined  in  Section  3.1  was  used  for  simulation  of 
one  cycle  of  consecutive  adsorption-desorption,  i.e.,  adsorption  during  the  first  part  of  the 
process,  followed  by  electrothermal  desorption  under  constant  voltage  supply  and  with  flow 
of  pure  innert.  This  process  can  be  defined  by  the  following  initial  conditions: 

t<  0:  Tg=Ts=Tgit=Tgot=Tgin=Tgo=Ta=Tp,  C=Clt=Cot=C,n=C*=O(q)=Cp=0 

and  the  following  time  changes  of  the  inlet  concentration,  voltage  and  inert  flow-rate: 

0<t<ts:  C;„=consteO,  U=  0  (6Qej=Q),  G=Gads  -  adsorption  step 

t>ts:  Cin= 0,  £/=conste0  (8Qep^0),  G=Ged~  electrothermal  desorption  step 
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Some  of  these  simulation  results  are  shown  in  Figures  17-20.  In  these  simulations  the 
switching  time  between  adsorption  and  desorption  ts  was  chosen  as  the  breakthrough  time,  i.e. 
the  time  after  which  the  outlet  concentration  reached  5%  of  the  inlet  concentration.  As  in  the 
real  TSA  process  [3],  the  flow-rate  during  desorption  was  much  lower  than  during  the 
adsorption  step.  The  presented  results  were  again  obtained  for  10  collocation  points. 

The  time  profiles  of  the  concentrations  in  both  phases  in  all  collocation  points,  as 
well  as  the  gas  phase  in  the  inlet  and  in  the  outlet  tube  are  shown  in  Figure  17,  and  the 
corresponding  temperatures  in  Figure  18.  The  influence  of  the  electrical  voltage  during  the 
desorption  step  is  shown  in  Figure  19,  and  the  influence  of  the  flow-rate  during  desorption  in 
Figure  20. 


Figure  17  Simulation  of  consecutive  adsorption-desorption  -  time  profiles  of  the  concentrations  in  10  collocation  points  in 
the  ads.  bed  (colored  lines)  and  in  the  inlet  and  outlet  tubes  ( — ),  for  GADS=0.\  mol/s,  GED= 0.03  mol/s  and  £7=10  V: 

a)  solid  phase;  b)  gas  phase 


(a)  (b) 


Figure  18  Simulation  of  consecutive  adsorption-desorption  -  time  profiles  of  the  temperatures  in  10  collocation  points  in  the 
ads.  bed  (colored  lines)  and  in  the  inlet  and  outlet  tubes  ( — ),  for  G^^O.l  mol/s,  GED= 0.03  mol/s  and  U=  10  V: 

a)  solid  phase;  b)  gas  phase 
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Figure  19  Simulation  of  consecutive  adsorption-desorption  -  influence  of  electric  voltage  on  the  concentrations  and 
temperatures  at  r=r2  (solid)  and  in  the  outlet  tube  (gas)  for  GADs=0A  mol/s  and  GED= 0.03  mol/s: 

a)  concentrations  i  b)  temperatures 
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Figure  20  Simulation  of  consecutive  adsorption-desorption  -  influence  of  the  flow-rate  during  desorption  on  the 
concentrations  and  temperatures  at  r=r2  (solid)  and  in  the  outlet  tube  (gas)  for  GADs=0A  mol/s  and  U=  10  V 

a)  concentrations  i  b)  temperatures 
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3.7.  Modeling  and  simulation  of  adsorption  with  reversed  flow 


In  Section  3.1  the  model  was  postulated  for 
the  annular  cartridge-type  adsorption  bed  with  radial 
flow  from  the  central  pipe  towards  the  outer  surface 
of  the  cartridge,  as  shown  in  Figure  3.  Nevertheless, 
in  the  system  shown  in  Figure  1,  during  adsorption 
in  some  of  the  cartridges,  the  gas  flow  was  in  the 
reverse  direction,  i.e.  from  the  periphery  of  the 
cartridge  towards  its  central  pipe  (the  upper 
cartridge  in  Figure  2)  [3],  For  a  single  cartridge 
configuration  considered  in  this  report,  in  this  case 
the  gas  is  introduced  through  the  outer  annular  tube 
and  it  leaves  the  system  through  the  central  tube. 

The  model  equations  and  some  simulation 
results  corresponding  to  this  case  are  given  in  this 
section.  The  schematic  representation  of  a 
differential  volume  element  of  the  cartridge,  used 
for  setting  up  the  model  equations,  is  given  in 
Figure  21.  The  assumptions  listed  at  the  beginning 
of  Sections  3  are  valid  for  this  case  as  well. 


Figure  21  A  differential  volume  element,  used 
for  modeling  of  adsorption  with  reversed  flow 


Model  equations: 

Compared  to  the  model  given  in  Section  3.1,  the  only  differences  are  in  the  model 
equations  for  the  central  and  the  outer  annular  tube,  and  in  the  sign  of  the  flow-terms  in  the 
material  and  energy  balances  for  the  gas  phase  in  the  adsorbent  bed.  Consecutively,  the  model 
is  obtained  as  the  following  set  of  PDEs  and  ODEs: 

Material  balance  for  the  solid  phase  within  the  adsorption  bed:  same  as  in  Section  3.1  - 
equation  (1) 


Material  balance  for  the  gas  phase  within  the  adsorbent  bed: 


dC 


dC 


s  dt  (2rnHeb)  dr 


+  (i„a)(C-C*)  =  ^|- 

r  or 


r 


v 


dC_ 

dr 


(46) 


Heat  balance  for  the  solid  phase  within  the  adsorption  bed :  same  as  in  Section  3.1  -  equation 

(3) 

Heat  balance  for  the  gas  phase  within  the  adsorbent  bed: 


8T 


dC  G(c  +  c  C)  dT 


p  (c  +c  C)  — —  +  p  c  T  =  — — — - — 

g  Pg  dt  Vg  pv  g  dt  2rnHEh 


dr 


-  + 


2rnHeh  dr 


r  dr 


,K 

dr 


(47) 


Material  balance  for  the  gas  phase  in  the  inlet  annular  tube: 

-pnHp =  GC„  -  GC„  +  kml (2VH )(1  - s„ XC’^  -  C„ )  +  2r,nH 8, (C|m;  -  C, )  (48) 
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Heat  balance  for  the  gas  phase  in  the  inlet  annular  tube: 


r=r\\ 


r=r2: 


(r32 -r22)nHPg  f_(cpg  +  cpvCit)Tgit]  =  G(cpg+cpvCin)Tgin-G(cpg+cpvCit)Tgit 
+hs2(2r2nH)(l-eb)(Tlr__r2 -Tgit)  +  h^nHs^T^ -Tgit)~  hwg(2r3nH)(Tgit  -Ta) 


(49) 


Material  balance  for  the  gas  phase  in  the  outlet  central  tube: 


r2nHpg 


dC^ 

dt 


=  GC 

i 


~GCot  +kml(2rpiH)(\-&b)(C* 


~Cot)  +  2 rpiH &b (C  r=n  -  Cot )  (50) 


Heat  balance  for  the  gas  phase  in  the  outlet  central  tube: 


fnHo  — r  (c  +  c  C  ,)T  ,1  =  G(c  +c  C\  )T ,  -G(c  +c  C  ,)T  , 

1  LV  PS  Pv  otJ  sot  J  v  pg  pv  \r=rp  g\r=rx  V  pv  old  got 


(51) 


+hsl(2rlizH)(\-eb)(T]r__n -T^  +  h^nHe^T^ -Tgot) 


Initial  conditions  -  same  as  in  Section  3.1  -  equation  (11): 
Boundary  conditions : 


Dt 


hs 


dr 


=  K(Ts-Tgot),  D, 


hg 


ST 


dr 


=  K,(T-T,),D.f  I 


dr 


=  C-C 


(52) 


-D, 


hs 


K 

dr 


=  hs2  (Ts~Tgit),  -D; 


hg 


dT 


dr 


git ' 


dC_ 

dr 


=  C-Cit  (53) 


These  model  equations  were  solved  in  an  analogous  way  as  the  ones  defined  in 
Section  3.1,  i.e.,  by  transforming  the  PDEs  into  approximate  sets  of  ODEs,  using  the  method 
of  orthogonal  collocations,  and  solving  the  resulting  set  of  nonlinear  27V+6  ODEs  (N  being  the 
number  of  internal  collocation  points),  numerically. 

The  model  presented  in  this  section  was  used  for  simulation  of  adsorption  with  flow 
from  the  periphery  towards  the  center  of  the  cartridge.  A  sample  of  the  simulated  time  profiles 
of  the  solid  and  gas  concentrations  in  10  collocation  points  and  in  the  inlet  and  outlet  tubes  is 
presented  in  Figure  22,  and  the  corresponding  temperatures  in  Figure  23.  Most  of  the  model 
parameters  were  the  same  as  in  Section  3.4  (listed  in  Table  1),  nevertheless,  some  of  the 
parameters  had  to  be  changed  in  order  to  avoid  convergence  problems.  Instead  the  values 
given  in  Table  1,  the  following  values  were  used:  &mi=A:m2=3xlO'7  mol/(cm2s),  hg2=2x  1 0‘2 
W/(cm2K ),  Ahs=1.2xl0'5  W/(cmK). 

The  simulated  concentration  and  temperature  profiles  obtained  for  two  different  gas 
flow  directions  have  similar  shapes,  i.e.  no  substantial  difference  of  the  two  processes  could 
be  observed.  A  somewhat  unusual  shape  of  the  curves  corresponding  to  the  collocation  points 
close  to  the  outer  surface  of  the  adsorption  bed  could  be  a  result  of  numerical  unstability,  but 
it  needs  further  analysis. 
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Figure  22  Simulation  of  adsorption  with  reversed  flow-  time  profiles  of  the  concentrations  in  10  collocation  points  in  the 
ads.  bed  (colored  lines)  and  in  the  inlet  and  outlet  tubes  ( — ),  for  G=0.06  mol/s  :  a)  solid  phase;  b)  gas  phase 


(a)  (b) 


Figure  23  Simulation  of  adsorption  with  reversed  flow-  time  profiles  of  the  temperatures  in  10  collocation  points  in  the  ads. 
bed  (colored  lines)  and  in  the  inlet  and  outlet  tubes  ( — ),  for  G= 0.06  mol/s  :  a)  solid  phase;  b)  gas  phase 
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4.  Mathematical  model  of  the  single-cartridge  adsorber  with  a  layered 
(nonhomogeneous)  adsorbent  bed 


The  model  presented  in  Section  3  was  based  on  the  assumption  that  the  adsorbent 
bed  can  be  treated  as  homogeneous.  Nevertheless,  in  reality,  the  ACFC  adsorbent  bed  under 
consideration  is  made  of  several  layers  of  tissue  made  of  activated  carbon.  Each  layer,  on  the 
other  hand,  is  a  complex  weaved  structure  made  of  yams,  each  of  them  consisting  of  a  great 
number  of  fine  carbon  fibers.  The  model  presented  in  this  section  takes  into  account  the 
existence  of  separate  layers  of  carbon  tissue  in  the  bed,  but  neglects  the  finer  structure,  i.e.,  it 
assumes  that  each  layer  itself  is  a  homogeneous,  porous  structure. 

A  simple  analysis,  based  on  ACFC  density  experimental  data,  justifies  using  such  a 
model.  E.g.,  for  American  Kynol  ACFC-5092-20  adsorbent,  the  porosity  of  the  material  itself 
(single  layer)  is  0.513  [10],  which  is  considerably  lower  than  the  average  porosity  of  the 
whole  adsorption  bed  which  is  -0.72  [3],  Based  on  these  values  and  the  data  about  the 
thickness  of  a  single  ACFC  layer  A/=0.55  mm  [10],  it  was  possible  to  calculate  the  gap 
between  two  layers:  A//=0.157  mm. 

A  schematic  representation  of  the 
layered  ACFC  radial-flow  adsorption  bed  is 
shown  in  Figure  24.  Only  the  model  for  gas 
inlet  through  the  central  tube  will  be 
considered.  The  model  for  the  reversed  flow 
can  be  obtained  in  an  analogous  way  as 
described  in  Section  3.7. 

In  postulating  this  model,  two  zones 
of  the  adsorption  bed  were  considered: 

ZONE  I  -  consisting  of  the  carbon  tissue 
layers,  and  ZONE  II  -  consisting  of  the  gaps 
between  the  layers.  Each  zone  consists  of  a 
number  of  sections,  (sections  of  ZONE  I 
and  ZONE  II  interchanging  regularly). 

For  an  adsorption  bed  made  of  n 
layers  of  tissue,  the  mathematical  model 
consists  of  the  material  and  energy  balances  for  the  solid  and  for  the  gas  phase  within  the 
ACFC  layer,  the  material  and  energy  balances  for  the  gas  phase  in  the  gaps  between  the  layers 
and  the  material  and  energy  balances  for  the  gas  phase  in  the  inlet  and  in  the  outlet  tube.  The 
assumptions  listed  in  Section  3  were  used  for  postulating  of  this  model  as  well.  In  addition,  it 
was  assumed  that  the  gas  in  the  space  between  two  layers  is  practically  mixed. 


Slika  24  A  schematic  representation  of  the  annular  - 
radial  flow  -  adsorbent  bed  made  of  layers  of  ACFC 


4.1.  Model  equations 

For  these  assumptions,  the  mathematical  model  of  an  adsorption  bed  made  of  n 
layers  consists  of  An  PDEs  (mass  and  heat  balances  for  ZONE  I)  and  2«+2  ODEs  (mass  and 
heat  balances  for  ZONE  II  and  the  inlet  and  outlet  tube). 
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For  sections  of  ZONE  I  -  ACFC  layers 

Adsorbate  balance  for  the  solid  phase: 


pc  =  (k  a)(Ck  -  O ,  *=1,3,5,. ..,2/i-l 


Adsorbate  balance  for  the  gas  phase: 


dt  2rnHs„  or  r  or  \  or 


Heat  balance  for  the  solid  phase: 


p.  K + <;„?» )  % = % + f  f r  -  r8<- )  (56> 

ot  ot  dV  ot  r  or  \  or  J 

*=1,3,5,.. .,2«-l 


//eat  balance  for  the  gas  phase: 


p  (c  +  c  CA—^  +  p  c  T , — i 

ps  pv  k>  dt  y  g  pv  gk  Qf 


dCjL=  +cpvQ)  GepvTgk  QCk 
dt  2rnHe„  dr  2nzHs„  dr 


Df  d  (  dT 


(ha)(Tsk-TA  +  ^- 


r  dr  [  dr 


*=1,3,5,. ..,2n-l 


For  sections  of  ZONE  II  -  Gaps  between  two  layers 


Adsorbate  balance: 


(rLi  -  rl )nHpg  ^  =  G  C,_,  I  -  GCt  +  knk  (2rknH){\  -  sc XC>  I  - ' Ct )  ■ + 


'*+1  dt 


kmk+1(2rk+lnH)(\sc)(C;+l  -Ck)  + 

'r~rk+ 1 

2rkTzHec(Ck-i\r=^  —  Ck)  +  2rk+lnHec(Ck+l\r=r^  —  Ck) 


Heat  balance: 


*=2,4,6,.  ,.2n-2 


(r*+i  "r*  +SvQ)^]  =  +cpv  -G(c/)g  +c,vQ)/;, 

+^(2r^//)(l-8c)(4_1^  - r#)  +  hgk (2rknHec)(Tgk_llr=rt  -Tgk)  + 

Kk+ 1  (2/*+1  -  SC)(4+1|^+1  -  Tgk )  +  *gi+1  (2rk+inH  sc  X^+I|r=/l+1  -  ^ ) 

*=2,4,6,.  ..,2n-2 
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For  the  inlet  and  outlet  tubes: 


Adsorbate  balance  for  the  gas  phase  in  the  central  inlet  tube: 


r'*np  ‘lt  =  GC“  ~GCil  L,  -c„)+2*i*ffMC,L,  ~c») 


(60) 


Heat  balance  for  the  gas  phase  in  the  central  inlet  tube: 


:2nHp  —\(c  +  c  CAT  ,1  =  G(c  +c  CAT.  -G(c  +c  CAT  , 

1  Vg  j  \_\^pg  ^pv^itd  git  J  ^  pg  pv  in  /  gin  ^  pg  ^pv^itd  git 


(61) 


+K\ (2r17t//’)(l  -  sc  ){Tsl,  -  Tgit  )  +  hgl  {2rpiHzc  )(Tgiir=r -TgU) 


g\\r=rx  git ) 


Adsorbate  balance  for  the  gas  phase  in  the  outlet  annular  tube: 


=GC. 


dt 


2/7—1  \r=r 


GCot  +  km2n  (2r2nnH)(\  -  sc )(C2n_1 , 


■CJ  + 


(62) 


2r2nizH&c{C2n_x 


r=r2n 


-cj 


Heat  balance  for  the  gas  phase  in  the  outlet  annular  tube: 


( rot  -  rL )*H\ Pg  Jt  [(cpg  +  cpvCot  )Tgot  ]  =  G(cpg  +  cpv  C2„_, )Tgln_lV=rin  -  G(cpg  +  cpvCot  )Tgot 


(63) 


+hs2S2r2nKH){\-&c){T2n_^ -TgJ  +  hg2X2r2nKH&AiTg2n_^ -TgJ- hwg(2rolnH)(Tgol -Ta) 


In  equations  (54-59)  the  subscript  k  associated  to  the  concentrations  and  temperatures 
denotes  the  section  in  the  adsorbent  bed,  starting  from  its  inner  surface.  The  odd  subscripts 
correspond  to  the  sections  of  ZONE  I,  while  the  even  ones  correspond  to  the  sections  of 
ZONE  II.  The  subscripts  associated  to  the  radial  coordinate  r  correspond  to  the  inner  and 
outer  diameters  of  the  sections:  e.g.  the  section  k  is  bounded  by  iy  and  r*+i,  r\  is  the  diameter 
of  the  central  inlet  tube  and  r2n+i  of  the  outlet  tube. 


Initial  and  boundary  conditions: 

Initial  conditions: 

t<0,  re[ri,r2n]:  Tg=Ts=Tgit=Tgot=Tgin=Tgo=Ta=Tp,  C=Cit=Cot=Cin=C*=®(q)=Cp  (64) 


Boundary  conditions: 
dCx 


r=rp.D 


dr 


,  dT. 

=  C,-C„,  Df  gl 


dr 


=  K^-TS„).  D, 


h.  ST„ 


dr 


=  KAT,t-TgJ  (65) 


r=rk+ 1,  for  A=l,3,5,...,2n-3: 


-D  aC* 


dr 


,  dT, 

—  c  _c  —  nh§ 

—  w  w+ 1 5  ~ 

or 


~  hgk+\  (Tgk  Tgk+i )  5  A 


hs  8Tsk 


dr 


=  66) 
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r=r2n: 


r=rk,  k=3,5,...,2n-l\ 

dC, 


D, 


dr 


8T 

=  C,-CM  D,*-* 
or 


=  hAT,-Tk_l),  D; 


hs  dTsh 


r=rk 


dr 


=  hATsk-T,  i)  (67) 


-D, 


dC2n-i 


dr 


=C  -C  -D 

^2n-l  '-'ot  ’  I 


dT 

hg  2_g2n-l 


dr 


—  hg2n(Tg2n-\  ^got)  ’ 


hs  dTsln_x 


dr 


=  K2niTs2n.l-Tgot)m 


4.2.  Numerical  solution  of  the  model  equations 


The  model  equations  presented 
in  this  section  were  again  solved 
numerically.  The  first  step  of  the 
numerical  solution  was  approximation 
of  the  PDEs,  using  the  method  of 
orthogonal  collocation.  This  method  was 
applied  to  each  section  of  ZONE  I.  For 
convenience,  the  radial  coordinate  was 
again  transformed  into  nondimensional 
form.  For  the  k-  th  section,  the 
nondimensional  radial  coordinate  was 
defined  in  the  following  way: 


i-0 


r2 


(69) 


Figure  25  A  schematic  representation  of  the  distribution 
of  the  collocation  points  in  the  first  3  sections 


For  each  section  of  ZONE  I,  N  internal  collocation  points  were  used,  plus  2  points  at  the 
section  boundaries.  (In  general  the  number  of  collocation  points  for  different  sections  could 
be  different,  but  in  the  analysis  presented  here,  the  same  number  was  used  for  all  sections.) 

As  an  illustration,  a  schematic  picture  of  the  distribution  of  the  collocation  points 
and  definition  of  the  variables  in  the  first  3  sections  (the  first  two  layers  of  tissue  and  the  gap 
between  them)  is  given  in  Figure  25. 

Application  of  the  orthogonal  collocation  method  on  PDEs  (54-57)  leads  to  the 
following  set  of  ODEs: 


dclk.i  Ka 


dt 


Pb 


(Q,-C;,),  1=0,.. .N+\,  k=l,3,5,...,2n-l 


N+ 1 


dC, 


k,i 


g2X,.A, 

j=0 


N+ 1 


ka 


D  Y  B,  C. 

k,i,j  k,j 


dt 


2TzHecA,pg(ukJA1+rk)  pg 
/=!,.. ..,N,  A^l,3,5,...,2«-1 


(Q,-c;,)+ 


7=0 


(Aiukj+rk)Pg 


(70) 


(71) 
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(72) 


dT 


sk,i 


U2 


-  +  - 


AH 


ads 


dq 


k,i 


dt  Pb(cps+cplqki)po0-  +  b(Tski-TR))H  (cps+cplqki)  dt 

(ha)(Tsk  i  -Tgk  i) 


+ 


N+\ 


+  - 


7=0 


cpiTsk,i  dqki 


P b (cPs  +  cpHk,i )(uk,Ai+rk)  p b (cps  +  cpHk,i )  ( cps  +  cpHk,i )  dt 

N,  k=\,3,5,...,2n-\ 


N+ 1 


7=0 


77+1 


GCpyTgkj  X  AjjCkj 
_ _ _ +o _ 

dt  2pg  (m/c,A,  +  ^  )nHscAj  pg  (cw  +  )2nHecA 
ha(Tsk  i  -Tgk  i) 


-  + 


7V+1 


-  +  - 


7=0 


j t  dCki 

Cpv  gKi  dt 


p g  (cPg + cpvck,i )  («*/ A/ + rk  )pg  (cPg  ~  cPvCk,i )  +  cpvQ,,- 

i=\,...,N,  k=  1,3,5,...,2«-1 


(73) 


In  equations  (70-73)  the  concentrations  and  temperatures  are  defined  by  two 
subscripts:  the  first  one  corresponds  to  the  section  and  the  second  one  to  the  collocation  point 
of  that  section. 

The  equations  corresponding  to  the  sections  of  ZONE  II  and  the  inlet  and  outlet 
tubes  have  also  be  rewritten: 


dCk  _  GCk_I  N+]  -GCk 
dt  (A  -rk)nHpg 


+ 


kmk(2rk)(l-ec)(C*  -Ct) 


(rL~rk)ps 


kmk+\  (2^+1  )(1  -  £c  )(C; 


k+ 1,0 


<++l  ~rk  )P* 


'^7)  |  2^£c(C/i._I  iV+i 


(^7+1  -  )Pe 


'  d'k )  +  2r4:+1sc(CJt+lj0  Cj.) 


<++,  -7.-  )p« 


£=2,4,...,2n-2 


(74) 


dT, 


G(c  +c  C, 

y  PS  PV  l 


pv^k-l,N+\ 


)T* 


GT 


gk 


hsk(2rk)(l-ec)(Ts 


sk-\,N+l 


~Tsk) 


dt  C rt+1  -  rk  )kH  pg  (cps  +  cCk )  (/£,  -  rk  )nH  pg  (rm  -rk)  pe  (cps  +  cCk ) 


-  + 


bgk(^-rk^c)^ek- 


gk-\,N+\ 


(n 


k+ 1  r*)p*(cw +<>£*) 


' d’gk )  +  Kk+i  (2rt+i )(1  £c)(^i:+i,o  'b'gk)  ^bgk+l(2rk+l£c)(Tgk+l0- 


Ak) 


(A 


rk  )P g(cPg+cPVCk) 


pv^kJ 


(£l 


-rk)p  g(cpg+cpvCk) 


(75) 


cpASk 


dC, 


(cPg+cpAk)  dt 

k=2,A,...,2n-2 


dCjL  = 
dt 


G(C j„  ~ Cjf )  |  (2/|Ti//)(l  sc)(C10  C,)  ^  2]]nHsc(C\  ()  Cit ) 

/f^T/p,,  p,,  K^HPg 


(76) 
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dT^  G(cre+cpvC,„)rg,  GT#  |  M^igXl-sJC T'^-Ty)  < 
dt  r2nHpg  (cpg  +  cpvCit )  /,2tc//  p?  7,2Ti//p?(c/)g +cn,C„) 

/7yl(2/;Ti//£c)(r  ,n-rir)  dc„ 

r2nH pg  (cpg  +  cpvCit )  (c/Jg  +  c/n,C„  )  dt 


d£ 0± 
dt 
(78) 


^(^2«-l,7V+l  Q)  ^m2n(^-r2n^^)(^  £c)(^2n-\,N+\  ^ ot )  '^r2n^^£ci^2n-\,N+\  ^ ot ) 


(r02-r22>//ps 


-  + 


(r02-r22>//ps 


+  - 


(rc2-r22>//ps 


dT 


got 


dt 


+  A^vC2n_1Af+1)7^2w_1>Ar+1 _ GTgn, _ ^  £c)(^2»-l,7V+l  ^gof)  | 

O'*  -  4  Pg  (CK  +  cpvCot )  irl  -  r2n  )nH pg  (r02  -  r2n  )nH pg  (cpg  +  cpvCot ) 


hg2n  (2r2 „;i7f sc )(rg2„_i  y+1  -  rgot )  (2roinH)(Tiot  -  Ta ) _ dCot 

(rot  -  rL )nH Pg  (cpg  +  cPvCot )  Oi  -  rl )nH pg  (cpg  +  cpvCot )  ( cpg  +  cpvCot )  dt 


(79) 


After  applying  the  method  of  orthogonal  collocation,  the  boundary  conditions  were 
obtained  in  the  algebraic  form: 


For  k=  1:  C,  0  = 


1  A+l 

c«+D.-jrL4».£i 

A/  J= 1 


1  +  0—4 


A, 


i,o,o 


T  = 

^gi,o 


1  7V+1 

"..V. 

I  j- 1 _  rr 

’  ^1.0  _ 


hJ.,+Dr — 
A 


i  7V+1 

“  rlA.,/. 

I  7=1 


si, 7 


h  _i_  J_  // 

"si  +  ^  a  A, 0,0 
A/ 


(80) 


For  ^=3,5,... ,2^-1:  Ck0  = 


J  7V+1 

Qr-i,jv+i  +  Dm  —  ^  4,o,Q, 

A/  7=1 


l-Z>—4 


k,  0,0 


T  = 

1sk$ 


1  7V+1  i  7V+1 

Vrf-l,*.l  +  A*’  T-  Z  A,.,/,. J  ‘A-U1.1  +  D!“  T-  Z  Vi 

A/  7=1 _  T  _ _ A/  7=1 

’  ^,0 


gk,j 


h  —  nhs  \  j 

%  ^  ^  ^,0,0 


h  —  nhs  _J_  a 

ngk  Ut  ^  ^£,0,0 


(81) 
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For k=l,39...92n-3:  C 


£,TV+1 


1 

Q+1,0  “Aw  “T- ^4;,TV+ljAj 

_ A7  7=0 _ 

1  +  Aw  ^k,N+\,N+\ 


i  TV  I  TV 

h  T  -  nhs  —  V  //  T  h  T  -  nh§  —  V  A  T 

nsk+l1gk+ 1,0  A  A^t  k,N+l,j1skJ  ngk+l1gk+\,0  A  A*i  kiN+XJ1  gkj 


A  j= 0 


l^,tv+i 


h  +  _J_  // 

A^£+l  ^  ^  ^  ^fc,TV+l,TV+l 


T 

’  ^2 


g£,TV+l 


_ Mg 

h  +  _J_  // 

'V+l  ^  ^  ^  ^,TV+l,TV+l 


(83) 


1  N 

r  _  n  ——V  A  C 

^ ot  A  A^l  ^n-l,N+l,j^2n-l,j 


Fork  2n-l:  C2w_1AA+1— ■ 


Ago 


1  +  Aw  ^  Aw-1,TV+1,TV+1 


L  g2n-l,N+l 


hg2n  Tgot  Dt§  ^  ^2n-\,N+\,jTg2n-\, 
_ A/  7=0 _ 

**2, +^^-4,-1 

A7 


1  ^ 

h  T  —  Dhs  ——  'S'  A  T 

ris2n  got  A  A^l  ^n-liN+XJ1  s2n-\,j 


ls2w-l,TV+l 


_ A  7=0 _ 

h  _i_  JL  j 

'*s2w  ^  ^  ^  ^w-^TV+^TV+l 


-(84) 


with: 


My  “ 


dl 


kj 


V  M  J 


Ku,  = 


uk,Ai+rk 


uk=uk,i 


4 


y2/4 


,  (In  j  j 


+  - 


dl, 


k,j 


\  du,  j 

\  k  Suk=uk, 


(85) 


4.3.  Simulation  of  adsorption  and  electrothermal  desorption 

The  derived  model  was  used  for  simulation  of  adsorption,  electroresistive  heating 
and  electrothermal  desorption.  The  simulation  was  performed  for  an  adsorption  bed  made  of  8 
layers  of  ACFC.  The  same  system  was  used  as  in  Section  3.  Owing  to  the  model  complexity 
and  the  large  number  of  equations  to  be  solved,  the  model  was  solved  for  3  collocation  points 
in  each  section  of  ZONE  1(1  internal  plus  2  border  points).  The  model  parameters  that  differ 
from  those  used  in  Section  3  are  listed  in  Table  2.  The  rest  have  the  values  given  in  Table  1. 

A  small  portion  of  the  simulation  results  is  given,  as  illustration. 

Figures  26  and  27  show  the  simulation  results  for  adsorption  on  an  initially  clean 
adsorption  bed.  The  time  profiles  of  the  concentrations  in  the  solid  phase  in  different  layers  of 
the  bed  (sections  of  ZONE  I)  are  presented  in  Fig.  26(a)  (red,  green  and  blue  lines  correspond 
to  different  collocation  points),  and  the  time  profdes  of  the  concentrations  in  the  gas  phase  in 
both  zones  (red,  green  and  blue  lines  correspond  to  different  collocation  points  in  the  sections 
of  ZONE  I,  while  the  dashed  cyan  lines  correspond  to  the  sections  of  ZONE  II),  as  well  as  in 
the  inlet  and  outlet  tubes  (dashed  black  lines),  are  given  in  Fig.  26(b).  The  corresponding 
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temperatures  are  shown  in  Figure  27  (the  solid  temperatures  in  Fig.  27(a)  and  the  gas 
temperatures  in  Fig.  27(b)). 

An  interesting  result  of  this  simulation  is  that  the  lines  corresponding  to  the 
concentrations  of  the  gas  phase  at  two  neighbouring  surfaces  in  two  ACFC  layers  practically 
overlap,  as  well  as  the  concentration  of  the  gas  between  those  two  layers.  A  similar  situation 
is  with  the  concentrations  in  the  solid  phase.  Nevertheless,  the  lines  corresponding  to  the 
temperatures  of  two  neighbouring  surfaces  are  distinct. 


Table  2.  Additional  model  parameters  used  for  simulation  of  adsorption,  electroresistive  heating  and 
electrothermal  desorption,  using  a  model  for  a  layered  adsorption  bed _ 


Variable 

Notation  and  value 

Adsorption  Heating/Desorption 

Outer  diameter  of  the  outlet  annular  tube 

rot= 3.55  cm 

Cloth  thickness 

Ay=0.055  cm 

Gap  between  two  layers  of  ACFC 

A//=0.0157  cm 

Cloth  porosity 

8c=0.513 

Mass  transfer  coefficient 

/c, ,,=0.0004  mol/(cm2s) 

Mass  transfer  coefficients  from  the  solid  in  zone 

I  to  the  gas  phase  in  the  inlet  and  outlet  annular 
tube 

km\=km\6=5xl0'5  mol/(cm2s) 

Mass  transfer  coefficients  from  the  solid  in  zone 

I  to  gas  in  zone  II 

km2=km3=. .  .=km  15=2x10"'  mol/(cm2s) 

Solid  to  gas  heat  transfer  coefficient  in  zone  I 

h= 4x  1  O'4  W/(cm2K)  h=6x  1  O'4  W/(cm2K) 

Heat  transfer  coefficients  from  solid  in  zone  I  to 
gas  in  zone  II 

/*s2=/*53=...=/*si5=lxl0-5  W/(cm2K) 

Heat  transfer  coefficients  from  solid  to  gas  in 
the  inlet  tube 

/?si=1.35xl0"5  W/(cm2K) 

Heat  transfer  coefficients  from  solid  to  gas  in 
the  outlet  annular  tube 

/?si6=8.7xl0_6W/(cm2K) 

Heat  transfer  coefficient  from  gas  in  zone  I  to 
gas  in  the  inlet  tube 

/?gi=8xl0"5W/(cm2K) 

Heat  transfer  coefficient  from  gas  in  zone  I  to 
gas  in  the  outlet  annular  tube 

hgi6= 2xlO'5W/(cm2K) 

Heat  transfer  coefficient  from  gas  in  zone  I  to 
gas  in  zone  II 

hS2=hgi=. .  .=/ig  15=6x1 0"5  W/(cm2K) 

Heat  diffusivity  of  the  solid  phase 

Dths=5  x  1  O'6  W/(cmK) 

Heat  diffusivity  of  the  gas  phase  in  zone  I 

Dthg=  lxl  O'6  W/(cmK) 

Radial  mass  transfer  dispersion  coefficient 

Dm= 9  x  1 0"6  mol/ (cms)  Dm= 5  x  1 0"5  mol/ (cms ) 

Adsorbent  bed  density 

pt=0.276  g/cm3 
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Ts  (K)  q  (mol/g) 


x  10 3  0)  x  104  (b) 


26  Simulation  of  adsorption  for  8  layers  of  ACFC  -  time  profiles  of  the  concentrations  in  all  sections  of  ZONE  I  and 
ZONE  II  and  in  the  inlet  and  outlet  tube  ( — ),  for  G=0.1  mol/s:  a)  Solid  phase;  b)  Gas  phase 


(a)  (b) 


Figure  27  Simulation  of  adsorption  for  8  layers  of  ACFC  -  time  profiles  of  the  temperatures  in  all  sections  of  ZONE  I  and 
ZONE  II  and  in  the  inlet  and  outlet  tube  ( — ),  for  G=0.1  mol/s:  a)  Solid  phase;  b)  Gas  phase 


The  simulated  time  profiles  of  the  temperatures  during  electroresistive  heating  are 
shown  in  Figure  28.  On  the  other  hand,  the  simulation  results  of  electrothermal  desorption  are 
presented  in  Figure  29  (the  concentration  profiles  in  different  sections  of  ZONE  I  and  ZONE 
II)  and  Figure  30  (the  temperature  profiles). 
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(a) 


(b) 


Figure  28  Simulation  of  electroresistive  heating  for  8  layers  of  ACFC  -  time  profiles  of  the  temperatures  in  all  sections  of 
ZONE  I  and  ZONE  II  and  in  the  inlet  and  outlet  tube  ( — ),  for  G=0.02  mol/s  and  U=  8V:  a)  Solid  phase;  b)  Gas  phase 


x  10'4  (a)  x  10'4  (b) 


Figure  29  Simulation  of  electrothermal  desorptions  for  8  layers  of  ACFC  -  time  profiles  of  the  concentrations  in  all  sections 
of  ZONE  I  and  ZONE  II  and  in  the  inlet  and  outlet  tube  ( — ),  for  G=0.02  mol/s  and  U=  8  V:  a)  Solid  phase;  b)  Gas  phase 


Figure  30  Simulation  of  electrothermal  desorptions  for  8  layers  of  ACFC  -  time  profiles  of  the  temperatures  in  all  sections  of 
ZONE  I  and  ZONE  II  and  in  the  inlet  and  outlet  tube  ( — ),  for  G=  0.02  mol/s  and  G=8V:  a)  Solid  phase;  b)  Gas  phase 
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5.  CONCLUSIONS 


This  report  summarizes  the  results  of  mathematical  modeling  of  the  central  unit  of  the 
TSA  system  with  electrothermal  desorption  -  the  adsorption  bed.  The  modeling  was 
performed  for  an  annular,  radial-flow  ACFC  adsorber,  used  in  Ref.  [3].  Two  models  of 
different  complexity  were  postulated  and  the  procedures  for  their  numerical  solution  were 
established.  The  models  were  used  for  simulation  of  adsorption,  electroresistive  heating  and 
electrothermal  desorption,  as  well  as  combined  adsortion-desorption.  The  equilibrium,  and 
part  of  the  transport  parameters  in  the  model  equations  were  calculated  using  correlations 
from  literature.  Nevertheless,  some  of  the  parameters  had  to  be  assumed,  as  correlations  for 
their  estimation  are  not  available. 

The  idea  of  this  modeling  is  to  incorporate  the  derived  models  for  the  adsorber  into  a 
mathematical  model,  which  would  describe  the  whole  TSA  system.  As  the  adsorber 
represents  the  main,  and,  from  the  modeling  point  of  view,  the  most  complex  part  of  the 
system,  the  development  of  its  models  is  crucial. 

One  of  the  important  issues,  that  have  to  be  resolved  in  the  further  work,  is  getting 
reliable  correlations  for  calculation  of  the  transport  parameters  in  an  adsorbent  bed  with 
fibrous  structure,  such  as  the  one  made  of  ACFC.  To  our  knowledge,  such  correlations  are  not 
available  in  the  literature.  In  order  to  define  them,  serious  experimental  investigation  is 
needed.  The  other  issue  is  the  necessity  for  experimental  verification  of  the  models,  which 
would  enable  their  comparison  and  critical  evaluation. 


Nomenclature 

a(m2/m3) 

-  Specific  surface  area 

b(K~l) 

-  Temperature  coefficient  of  the  bed  electrical  resistivity 

C(mol/mol) 

-  Adsorbate  concentration  in  the  gas  phase 

C*(  mol/mol) 

-  Adsorbate  concentration  in  the  gas  phase  in  equilibrium  with  the  solid  phase 

cpg(J/mol/K) 

-  Specific  heat  capacity  of  the  inert  gas 

cpl(J/mol/K) 

-  Specific  heat  capacity  of  liquid  adsorbate 

cps(J/g/K) 

-  Heat  capacity  of  the  solid  phase 

cpv(J/mol/K) 

-  Specific  heat  capacity  of  gaseous  adsorbate 

Dab(  mol/cm/s) 

-  Molecular  diffusivity  of  the  gas  phase 

Z)m  (mol/cm/s) 

-  Radial  mass  transfer  dispersion  coefficient 

A>(cm) 

-  Fiber  diameter 

Ahg(W/K/cm) 

-  Heat  diffusivity  of  the  gas  phase 

Ahs(W/K/cm) 

-  Heat  diffusivity  of  the  solid  phase 

AJ/mol) 

-  Adsorption  energy  of  the  adsorbate  (D-R  eq.) 
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g(cm/s2) 

-Gravitation  constant 

G(mol/s) 

-  Flow  rate  of  the  inert  gas 

//(cm) 

-  Bed  axial  dimension 

/?(J/cm2/K) 

-  Solid  to  gas  heat  transfer  coefficient 

/?wg(J/cm2/K) 

-  Gas  to  ambient  heat  transfer  coefficient  (heat  losses) 

/(W/cm/K) 

-  Gas  phase  thermal  conductivity 

4(J/K) 

-  Boltzmann’s  constant 

/m(mol/cm2/s) 

-  Mass  transfer  coefficient 

M4(g/mol) 

-  Adsorbate  molar  mass 

Ms(g/mol) 

-  Inert  gas  molar  mass 

N 

-  Number  of  internal  collocation  points  in  the  radial  direction 

Nu 

-  Nusselt  number 

p(kPa) 

-  Gas  pressure 

Pa 

-  Adsorbate  partial  pressure 

/?cA(kPa) 

-  Critical  pressure  of  the  adsorbate 

/?cB(kPa) 

-  Critical  pressure  of  the  inert  gas 

//(kPa) 

-  Adsorbate  saturation  vapor  pressure 

Pe 

-  Peclet  number 

Pr 

-  Prandtl  number 

q(  mol/g) 

-  Adsorbate  concentration  in  the  solid  phase 

2el(J/s) 

-  Joule’s  heat 

Ra 

-  Rayleigh  number 

Re 

-  Reynolds  number 

/?g(J/mol/K) 

-  Gas  constant 

ri(cm) 

-  Inner  diameter  of  the  cartridge 

r(cm) 

-  Bed  radial  dimension 

/?el(^) 

-  Electric  resistance 

Sc 

-  Schmidt  number 

Sh 

-  Sherwood  number 

UK) 

-  Ambient  temperature 

/ca(K) 

-  Critical  temperature  of  the  adsorbate 

TcbCK) 

-  Critical  temperature  of  the  inert  gas 

T’g(K) 

-  Gas  phase  temperature 

?r(K) 

-  Referent  temperature 

n(K) 

-  Solid  phase  temperature 

C/(V) 

-  Electric  potential  difference 

U 

-  Non-dimensional  radial  column  coordinate 
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V 

VP  A,  VPB,  VP \ 
ff(cm3/g) 

W0(  cm3/g) 
A^ads(J/g) 

Greek  letters: 

a(cm2/s) 

P(l/K) 

A(cm) 

£b 

£c 

p  (Pa/s) 

v(cm2/s) 

p(Qcm) 

po(ficm) 

pb(g/cm3) 

pg(mol/cm3) 

pA(g/cm3) 

Subscripts: 

b 

c 

g 

in 

o 

P 

s 

it 

ot 

I 

II 


-  Superficial  gas  velocity 

,  VPd  -  Wagner  constants  for  the  adsorbate 

-  Adsorbate  volume  per  lg  of  the  adsorbent  mass 

-  Total  volume  of  micropores  (D-R  equation) 

-  Heat  of  adsorption 


-  Thermal  diffusivity  of  the  gas  phase 

-  Reciprocal  of  the  film  temperature 

-  Bed  thickness 

-  Bed  porosity 

-  Cloth  porosity 

-  Dynamic  viscosity 

-  Kinematic  viscosity  of  the  gas  phase 

-  Electric  resistivity 

-  Electric  resistivity  at  referent  temperature  TR 

-  Adsorbent  bed  density 

-  Density  of  the  inert  gas 

-  Adsorbate  density 


-  bed 

-  cloth 

-  gas 

-  inlet 

-  outlet 

-  previous  (initial) 

-  solid  phase 

-  inlet  tube 

-  outlet  annular  tube 

-  ZONE  I 

-  ZONE  II 
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